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The transfer matrix formalism is implemented in the form of the multiple collision technique 
to account for dissipative transmission processes by using complex potentials in several models of 
atomic chains. The absorption term is rigorously treated to recover unitarity for the non-hermitian 
hamiltonians. In contrast to other models of parametrized scatterers we assemble explicit potentials 
profiles in the form of delta arrays, Poschl- Teller holes and complex Scarf potentials. The techniques 
developed provide analytical expressions for the scattering and absorption probabilities of arbitrarily 
long wires. The approach presented is suitable for modelling molecular aggregate potentials and 
also supports new models of continuous disordered systems. The results obtained also suggest the 
possibility of using these complex potentials within disordered wires to study the loss of coherence 
in the electronic localization regime due to phase-breaking inelastic processes. 
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I. INTRODUCTION 

The inelastic scattering processes occurring in meso- 
scopic samples as a consequence of a finite non-zero tem- 
perature can noticeably change the coherent transport 
fingerprints of these structures. The worsening of elec- 
tronic transmission due to such effects is expected but 
in some situations the competition between the phase- 
breaking mechanisms and the quantum coherent inter- 
ferences can improve conductance in certain energetic 
regimes. This is the case, for example, of disordered 
structures. This fact has attracted much attention in 
the study and modelling of dissipative transport in one- 
dimensional structures. Interest is also prompted by 
experiments currently being carried out on real atomic 
chains 0. 

A model of parametrized scatterers coupled through 
additional side channels to electron reservoirs incorpo- 
rating inelastic events was initially proposed by Biittiker 
[^], and much work has been done along this line 
On the other hand, inelastic processes can be modelled 
by small absorptions which in turn can be described by 
extending the nature of the quantum potentials to the 
complex domain. The main purpose of this work is to in- 
clude absorptive processes by performing these complex 
extensions on previous quantum wire models developed 
by the authors 3 and also on other atomic potentials. 

The use of complex site energies and frequencies has 
already been considered in the study of electronic con- 
ductivity through one-dimensional chains 0, Q , but non- 
hermitian hamiltonians have also been used to account 
for a large variety of phenomena, ranging from wave 
transport in absorbing media 7|, violation of the single 
parameter scaling in one-dimensional absorbing systems 
[Sj , appearance of exceptional points in scattering theory 
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and quantum chaology description of vortex de- 
localization in superconductors with a transverse Meiss- 
ner effect [Tll | and more phenomenologically with nuclear 
optical potentials. Special mention is required for the 
framework of T^T-symmetry [l^ . where it is possible to 
consider periodic wires under complex potentials show- 
ing real band spectra [TsL IT^ . There is nothing wrong 
in principle with the use of non-hermitian hamiltonians 
as long as their properties are controlled by a sufficient 
knowledge of the full spectrum. Indeed, renormalization 
group calculations have been carried out giving rise to 
imaginary couplings as a result of quantum dressing of 
the classical real potentials 

An interesting modern review on absorption in quan- 
tum mechanics has appeared recently and we address 
the interested reader to this publication and references 
therein. 

The paper is organized as follows. In Section ^ we 
briefly review the multiple collision technique based on 
the transfer matrix method, and in Section IlIII we show 
how unitarity can be easily restored in the presence of 
absorption and how the general unitarity condition can 
be generalized accordingly. We then turn our attention 
to arrays of delta potentials and calculate and draw the 
scattering and absorption probabilities in Section llVI In 
Section the Poschl- Teller potential is used to build 
atomic chains, and its complex extension the complex 
Scarf potential is fully developed in Section IVII The 
analytical scattering probabilities are shown for a vari- 
ety of composite potential profiles and the effect of the 
imaginary parts on the transmission is analyzed. The 
calculations concerning exact wave functions and analyt- 
ical conditions are offered in two appendices. The paper 
ends with several concluding remarks. 



II. THE MULTIPLE COLLISION TECHNIQUE 

The time-independent scattering process in one dimen- 
sion can be described using the well known continuous 
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transfer matrix method [l7l| . 
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where Al,Bl{Ar,Br) mean the amphtudes of the 
asymptotic traveUing plane waves g^^x^^-ikx ^^^^ j^j-^ 
(right) side of the potential. Whatever the nature of the 
potential is, real or complex, the transmission matrix al- 
ways verifies det M = 1 as a consequence of the constant 
Wronskian of the solutions of the Schrodinger equation. 
The transmission and reflection amplitudes then read, 
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where the superscripts L, R stand for left and right inci- 
dence. The insensitivity of the transmission amplitude to 
the incidence direction is a universal property that holds 
for all kind of potentials. However, the reflectivity, al- 
though symmetric for real potentials, changes with the 
incidence side for a complex one unless it is symmetric 
[isj l . The effect of a composition of n different potentials 
can then be considered as the product of their transmis- 
sion matrices, 



M„ . . .M2M1. 



(3) 



The transmission matrix formalism is an important tool 
for the numerical treatment of different problems. An 
intuitive and general interpretation of the composition 
procedure can be given in the following form. Consider 
two potentials Vi (x) , V2 (x) characterized by the scatter- 
ing amplitudes ti, , r^, t2, ^2 , and joined at a certain 
point. Then, the scattering amplitudes of the composite 
potential can be obtained by considering the coherence 
sum of all the multiple reflection processes that might 
occur at the connection region. 
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Replacing the scattering amplitudes with the elements 
of the corresponding transmission matrices Mi , M2 , one 
can trivially check that in fact these last formulae are 
the equations of the product M2M1. Thus, the compo- 
sition rules given by Q are not restricted to the conver- 
gence interval of the series X]^o('^2'^?)"- They provide 
an explicit relation of the global scattering amplitudes 
in terms of the individual former ones and can be easily 
used recurrently for numerical purposes. This composi- 
tion technique was first derived for a potential barrier 
pro] and has been used for designing absorbing potentials 



III. THE SCHRODINGER EQUATION FOR A 
COMPLEX POTENTIAL 



Let us consider a one-dimensional complex potential 
of finite support Vix) = Vr{x) + iV^{x) {V{±oo) = 0). 
For the stationary scattering states, the density of the 
current flux is proportional to the imaginary part of the 
potential 

^ = V,(.)|vl/(x)p (5) 
ax n 

where J{x) is defined as. 

Therefore, in the presence of a non- vanishing Vi{x) the 
unitarity relation regarding the transmission and reflec- 
tion probabilities T{E) + R{E) = 1 is no longer valid. 
One can still recover a pseudounitarity relation by defin- 
ing a quantity that accounts for the loss of flux in the 
scattering process. Dealing with the asymptotic state 



-ikx 



, *^(+oo) 



t{k)e 



ikx 



can write the asymptotic values of the flux as, 
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This latter equation remains the same for the right inci- 
dence case (with R^(k)) when the asymptotic state takes 
the form *^(-oo) = t{k)e-'''^ , (+00) 
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Using eq. lO the flux term reads. 
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and it is usually understood as the probability of absorp- 
tion But A{k) must be a positive defined quantity 
in order to be strictly considered as a probability and this 
is not ensured by the definition (unless Vi{x) < Vx). 
The sign of A{k) depends on both the changes in sign of 
the imaginary part of the potential and the spatial distri- 
bution of the state. Although a negative value for A{k) 
could be viewed as emission (because it means a gain in 
the flux current) it also leads the transmitivity and the re- 
flectivity to attain anomalous values T{k) > 1, R{k) > 1, 
which are difhcult to interpret. Let us also note that the 
integral representation of the absorption term is useless 
for practical purposes because to build the correct ex- 
pression of the state \I'^'^(a;) one needs to impose the 
given asymptotic forms to the general solution of the 
Schrodinger equation, therefore obtaining the scattering 
amplitudes, so one cannot calculate the absorption prob- 
ability without knowing R{k) and T{k). 
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IV. SCATTERING OF A CHAIN OF DELTA 
POTENTIALS 



matrix for a delta potential preceded by a zero potential 
zone of length a reads, 



Let us consider a potential constituted by a finite ar- 
ray of Dirac delta distributions, each one with its own 
coupling at and equally spaced at a distance a. This 
is probably the simplest one-dimensional model imagin- 
able, but in spite of its apparently simplicity it supports 
an unexpected physical richness. It has been successfully 
used to model band structure in a periodic quantum wire 
PI and has proved its usefulness when considering uncor- 
related and correlated disorder structures [3,|2l|, showing 
interesting effects such as the fractality of the density of 
states and the different localization regimes for the elec- 
tronic states. 

The global potential will be characterized by the ar- 
ranged sequence of the parameters (a/oj), where — 
means the "effective range" of the ith delta, in the 
order they appear from left to right. The transmission 
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Considering a chain of N different deltas and applying 
the composition rules to this type of matrices one finds 
that it is possible to write a closed expression for the 
scattering amplitudes. They are given by. 
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where for each j the means we are summing over 
the (^) combinations of size j from the set {1,2,..., N} 

a = {(Ti, (T2, . . . , CTj} with (Ti < (72 < . . . < (Tj. The 
amplitude up to a phase is obtained from for the re- 
verse chain. These latter formulae resemble the equations 
for the band structure and eigenenergies of the closed 
system [4|. In spite of their formidable aspect, Eas. (|13|) 
are easy to program for sequential calculations, provid- 
ing the transmitivity and reflectivity of the system with 
exact analytical expressions. 

Let us now incorporate the dissipative processes that 
are always present in real wires, causing energy losses. 
We have modelled that effect by including an imagi- 
nary part in the potential. In this case the natural 
complex extension of our system consists in promoting 
the delta couplings from real to complex, thus writing 
(a/flj) — rj — isj. We also take Sj > for all j in order 
to avoid anomalous scattering. The effect of including 
complex couplings on the spectrum of an infinite peri- 
odic delta chain has recently been studied in detail [T3| . 
Let us see what happens in a chain with open bound- 
aries. In Fig|21 the usual scattering diagram is shown 
for a short periodic chain with real potentials. Includ- 



ing a small imaginary part in the couplings we see how 
the transmission pattern is altered with a non-negligible 
absorption that peaks at the incoming band edges while 
the reflectivity is not noticeably changed. This tendency 
of the absorption term also appears when several species 
are included in the periodic array, and its pattern does 
not change much if only some of the couplings are com- 
plexified. 

When the array presents no ordering at all, the graph 
is quite unpredictable and different configurations can be 
obtained. In Fig^a) a peaky spectrum with very sharp 
absorption resonances is shown. The scattering process 
in this case is strongly dependent on the direction inci- 
dence, as can be seen. On the other hand, smoother dia- 
grams are also possible in which the effect of the complex 
potential manifests through an almost constant absorp- 
tion background and a small change depending upon the 
colliding side, like the one in Fig^b). 

This naive potential, apart from being exactly solvable, 
is powerful enough to account for very different physical 
schemes, which makes it a very useful bench-proof struc- 
ture. 
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FIG. 1: (Color online) Scattering process for disordered arrays of 15 deltas with complex couplings. The se- 
quences of the real parts of the characteristic parameters are: (a) Re(a/a-,): 3,1,2,0.5,3,2,1,3,0.5,4,5,1,2,2,3 (b) 
Re(a/aj): — 1, —4, —3, —1, —2, —3, —4, —1, —2, —3, —1, —4, —4, —2, —3 . The imaginary part of each coupling has been chosen as 
lm{a/aj) — — 0.01|Re(a/aj)|. The arrows in the legends mark the direction of incidence. 




FIG. 2: (Color online) Scattering and absorption probabilities 
for one-species delta chains with length A'' = 15 and param- 
eters (o/ai) — 1.0 (dashed lines) and (a/ai) — 1.0 — 0.015j 
(solid lines). 



V. ATOMIC QUANTUM WELLS 

Let us go one step further and consider a potential that 
resembles the profile of an atomic quantum well with an- 
alytical solutions, the well-known Poschl- Teller potential 
hole. It reads, 
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and it is shown in Fig|3| 

The probability of asymptotic transmission is well 
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FIG. 3: Poschl- Teller potential (arbitrary units) 
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One characteristic feature of the Poschl- Teller hole is that 
it behaves as an absolute transparent potential for inte- 
ger values of A, as can be seen from Eg. ljlSI) . From the 
wave functions one can also obtain the asymptotic trans- 
mission matrix, 
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where 7 = 2argr(iA;/a) - arg{r(A + ik/a)T{l - X + 
ik/a)}. Let us try to build a chain with these atomic 
units. In order to do so, one has to include a sensible cut- 
off in the potentials to ensure first that the wave function 
takes a proper form at the junction regions and second 
that the resulting potential hole can still be described 
by a handy transmission matrix, so that Eqs.lQJ can be 
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FIG. 4: (Color online) Transmission through a double Poschl- 
Teller hole with parameters a\ = =2 (a::"^ units), Ai — 
A2 = 2.4, di = d2 = 5 (a; units). The dashed line corresponds 
to a single potential hole. The inset shows the composite 
potential profile (arbitrary units). 
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FIG. 5: (Color online) Transmission patterns for two sym- 
metric composite potentials of ten units each. Their pro- 
files are shown in the insets (arbitrary units). Parameters 
for the first five potentials of the sequences: (solid line) 
a = 1, 1,0.5, 1, l(a;~^ units), A = 1.66,2.19,5.01,2.33,2.16, 
d = 4,4:, 7, 4, 4{x units) and (dashed line) q = 2, 2, 1, 1, l(x~^ 
units), A = 1.66, 1.66, 2.03, 2.03, 2.03, d = 2, 2, 4, 4, 4{x units). 



applied easily. The cut-ofF will be placed at a distance d 
from the center of the potential (Fig|31). The wave func- 
tion in the interval [—d, d] is ^'2(2:) = A2e(x) + B2o{x) 
where e{x),o{x) are the even and odd solutions respec- 
tively of the Schrodinger equation. Outside that interval 
the wave function is assumed to be a superposition of 
the free particle solutions, regions 1 and 3 in FigO The 
connection equations at the cut-off points lead to a rela- 
tionship between the amplitudes of the wave function in 
sectors 1 and 3 in terms of the values of e{x),o{x) and 
their spatial derivatives at ±d. Therefore, the distance d 
must be such that the asymptotic form of the solutions 
of the Schrodinger equation can be used at that point in 
order to ensure a sensible transition to the free particle 
state and to obtain a transmission matrix as simple as 
possible. The solutions e{x), o{x) as well as their asymp- 
totic forms are found in Reference nevertheless they 
are also reproduced in Appendix IXI 

After some algebra one finds the transmission matrix 
for the cut-off version of the potential hole, 
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ip 



-le 



(17) 

The matrix is the same as for the asymptotic case in 
Eq. (|16() plus an extra phase term in the diagonal elements 
that accounts for the distance 2d during which the parti- 
cle feels the effect of the potential. These phases are the 
key quantities for the composition procedure since they 
will be responsible for the interference processes that pro- 
duce the transmission patterns. Due to the rapid decay 
of the Poschl- Teller potential, the distance d admits very 
reasonable values. In fact, we have seen that for a sensi- 
bly wide range of the parameters a S [0.1,3], A S [1,5] 
one can take as a minimum value for the cut-off distance 
do ~ 2dif2 — 3.5/q!, where (ii/2 is the half-width (Fig|2Jl. 
Taking d > do the connection procedure works really 
well, as we have checked in all cases by comparing the 
analytical composition technique versus a high precision 
numerical integration of the Schrodinger equation, ob- 
taining an excellent degree of agreement. 

Composing two potential holes and applying Eqs.Q) 
together with Eq. H17() one finds for the transmission prob- 
ability, using the previously defined quantities p,j 



T2PT{k) = — , (18) 

pIpI + (1 +p?)(1 +pi) - 2pip2y^l+pWl+pi cos(7i + 72 + 2k{di +d2)) 



r 



which is a handy expression that can hardly be obtained not been made before. Ea. l|18|l clearly shows the inter- 
by trying to solve the Schrodinger equation for the dou- fcrcncc effect depending on the distance di + d2 between 
blc potential hole. To our knowledge this calculation has the centers of the holes. An example of transmission is 
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FIG. 6: Real and imaginary parts of the complex Scarf po- 
tential (arbitrary units). 

shown in Fig^ 

The composition procedure can be applied with a small 
number of atoms to study the transmitivity of different 
potential profiles resembling molecular structures such as 
those in Fig|31 The transmission matrix H17|) can also be 
used to consider a continuous disordered model in the 
form of a large chain of these potential holes with ran- 
dom parameters. So far, in the literature only two kinds 
of potentials have been used to build continuous disor- 
dered models, namely the Dirac delta potential and the 
square well(barrier) , due to their well known and easy 
to manipulate transmission matrices. We recall the fact 
that handy transmission matrices can be obtained for 
other potential profiles using reasonable approximations, 
such as the one described here. 

The next step for our purpose is to consider dissipation 
in these one-dimensional composite potentials. 



VI. DISSIPATIVE ATOMIC QUANTUM 
WELLS /BARRIERS 

We shall consider the extension of the Poschl- Teller 
potential given by the complexified Scarf potential, 

V{x) = ^ +iV2 J — -] (19) 

2rn Vcosh {ax) cosh (ax) J 

with Vi,V2 G M. It is a proper complex extension for two 
reasons: it admits analytical solutions ^23] and its imagi- 
nary part is somehow proportional to the derivative of the 
real potential. This latter criterion has been considered 
in nuclear optical potentials to choose adequate complex 
extensions. It seems reasonable to measure the strength 
of the dissipation processes in terms of the "density" of 
the real interaction and therefore writing an imaginary 
potential that is proportional to the spatial derivative of 
the real one. The potential profile is shown in FigEl 

The Scarf potential has been extensively considered in 
the literature, mainly dealing with its discrete spectrum, 
either in its real and complex forms, from the point of 
view of SUSY Quantum Mechanics |2J] or focusing on its 
T'T-symmetric form |25|. 

First, a detailed mathematical analysis of the poten- 
tial, regarding its scattering properties, must be made to 
discuss some new features and some assertions that have 
been made. 

The left scattering amplitudes of the real Scarf po- 
tential have been obtained in terms of complex Gamma 
functions |2^ . Recently, a considerable simplification has 
been pointed out by Ahmed '25]. In fact, the asymptotic 
transmitivity and reflectivity for the complex Scarf can 
be written as. 



T{k) = 



sinh^(27rfc/a) 



R^{k) 



sinh (27rfc/a) -I- 2cosh(27rfc/Q;) cosh(7r(7+) cosh(7r5_) -I- cosh (7rg+) + cosh (ng^) 

cosh^(7r5+)e-2'^'=/" + cosh^(7r.9_)e2'^'=/" + 2cosh(7r5+) cosh(7r5_) 
sinh^(27rfc/Q!) + 2 cosh(27rfc/a) cosh(7rg+) cosh(7rg_) + cosh^(7r5+) + cosh^(7r5_) 



(20) 
(21) 



where g± = \/Vi ±¥2 — 1/4: and R^{k) is recovered from 
R^{k) by interchanging and g_ (which is equivalent 
to substituting V2 —V2 and therefore changing the 
direction of incidence) . These expressions derive from the 
asymptotic transmission matrix, which is obtained here 
using the asymptotic form of the Schrodinger equation 
solutions (Appendix IB|) . 



Ms^arfik)^^ ^_ (22) 



where 

_ cosh(7r5+)e'''=/" +cosh(7r5_)e"^'=/" 
^ ^ sinh(27rfc/Q;) ^ ' 

s^sig+^g^) (24) 

and ip = 2arg{r(ifc/a)r(l/2 + ik/a)} - arg{r(c + 
ik/a)T{b + ik/a)T{l - c + ik/a)T{l -b + ik/a)} with 
the definitions 

c=l-^{9+-9-) ; b=^-^{g++g^). (25) 
It immediately follows from the transmission matrix that 
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\V2\ 



V2I emissive 



(-1^1,1^^21) 



-0.5 [0, 0.5] U [1.5, 5.5] U [11.5, 19.5] , 
-1 [0, 5] U [11, 19] U [29, 41]... 

-2.4 [0, 0.4] U [2.4, 3.6] U [9.6, 17.6] . 
-3 [0, 1] U {3} U [9, 17] . . . 

-4 [0, 4] U [8, 16] U [26, 38] . . . 
-5 [0, 1] U [3, 5] U [7, 15] . . . 



[0,0.5] 


[0,0.4] U [2.4,2.569] 

[0,1] 
[3.606, 4] 
[3,4.123] U [4.472,5] 



TABLE I: Ranges of V2 compatible with the condition T{k) < 
1 V fc for the complex Scarf potential for certain negative 
values of Vi . The last column includes the intervals providing 
physical scattering from the emissive side of the potential. 



the absorption probabilities read. 



1 



ss — 
l + ss' 



(26) 



Unlike the complex delta potentials example this po- 
tential has some drawbacks that must be carefully solved. 
Its imaginary part is non-negative defined in its domain, 
which might cause anomalous scattering. Only some val- 
ues of V2 will be physically acceptable. To ensure that 
T{k) < 1 Vfc, it is clear from Ea. H20l) that the necessary 
and sufficient condition is cosh(7r(7+) cosh(7rg_) > 0. The 
functions 5+, g- can be real or pure imaginary depend- 
ing on the values of Vi and V2. A detailed analysis of 
the conditions for physical transmission is presented in 
Appendix ^ As a summary, let us say that for Vi > Q 
(barriers) , the evaluation of the condition translates into, 

m e[o,v-i]u 

[2n{2n - 1) + Vi,2n{2n+l) + Vi] n e Z+ . 

(27) 

For Vi < (wells) the situation becomes more compli- 
cated and the result can only be expressed through sev- 
eral inequalities, each one adding a certain allowed range 
for V2 (see Appendix IB|) . As an example, in Table we 
show the compatible ranges of V2 for a few negative values 
of Vi . One can trivially check the compatibility of the in- 
tervals presented for V2 with the condition T{k) < 1 V fc 
by plotting Ea. l(^ . In a two dimensional plot of IV2I vs. 
Vi , the physical ranges for the transmission distribute as 
alternating fringes and a funny chessboard like pattern 
(Fig®. 

One feature to emphasize according to the conditions 
given for acceptable transmission is the fact that the 
number of permitted intervals for V2 is infinite for any Vi , 
either positive (barrier) or negative (well), and therefore 
there is no mathematical upper bound on | V2 1( | ^2'^''***'^"' |) 
above which the transmission probability always becomes 
unp hysical, contrary to what has been reported recently 
|25|. From a physical viewpoint of course, a sensible lim- 
itation must also be imposed on V2, usually IV2I ^ \Vi\. 

Let us see now what happens with the reflectivity. We 
assert that for the values of Vi and V2 preserving a phys- 
ical transmission, one of the reflectivities of the system 



(1,1) 
[2,0] 
[6,0] 
[12, 0] 



(3,3) 
(4,2) 
(9,3) 
(16,4) 



(6,6) 
(7,5) 
(13,7) 
(21,9) 



(10, 10) 
(11,9) 
(18, 12) 
(27, 15) 



(15,15) 
(16, 14) 
(24, 18) 
(34, 22) 



TABLE IL Some correlated values of Vi,V2 producing the 
fully resonant behaviour of the complex Scarf potential. The 
particular Poschl- Teller resonances are in square brackets. 



remains physical (i.e. R{k) < 1 — T{k) Vfc), left or 
right, depending on the particular values of Vi,V2 (or 
equivalently, one of the absorptions takes positive values 
for all fc). The statement is easy to prove from Eas. H26() 
and more specifically reads: when T{k) < 1 Vfc (i.e. 
cosh(7rg+) cosh(7r,g_) > 0), then 



I cosh(7r5_)| < I cosh(7r(;+)| ^ A^(fc) > Vfc 
|cosh(7r5_)| > |cosh(7r.g+)| ^ A'^ik) > Vfc. 



(28) 



Considering Vi > and V2 > it is not hard to see that 
the first of the above inequalities always holds. There- 
fore, in the case of a potential barrier the scattering is 
always physical from the absorptive side (trough of the 
imaginary part), as has already been stressed [23. More 
interesting is the fact that this conclusion cannot be ex- 
tended to the case Vi < (well). In this case the phys- 
ical scattering sometimes occurs from the emissive side 
(peak of the imaginary part), producing smaller absorp- 
tion terms. In Table HI a few examples of V2 intervals 
providing physical scattering from the emissive side for 
some potential wells are shown. 

Another interesting feature that must be observed is 
that there exists a set of correlated values of Vi,V2 for 
which the complex Scarf potential behaves as fully trans- 
parent. The condition for this to happen is from Ea. l|2U|) 
cosh(7rg-|_) = cosh(7r(7_) = 0. Thus, the main require- 
ment is that g+ , g_ must be pure imaginary, yielding in 
this case the transparency equations. 



(29) 



whose solution is, 

Vi = i-i[(2m+l)2 + (2n+l)2] 

V2 = ^ [(2m + 1)2 - (2n + 1)2] m,neZ. (30b) 
8 

It is worth noting that the transparencies only appear 
for potential wells (Vi < 0). Considering the particular 
case n ~ m one recovers the Poschl- Teller resonances 
(A e Z). In Tableinithc first values of Eqs. 1(201) are listed 
explicitly. 

The absorption obviously vanishes for all fc when con- 
sidering these special resonant values of the potential am- 
plitudes. Surprisingly, there also exists another set of 




(30a) 
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7: (Color online) Characteristic scattering patterns for 
complex Scarf well a = \{x~^ units), Vi = — 0.5, V2 = 
and (b) a complex Scarf barrier a = \{x~^ units), Vi = 
= 0.1. 



non-trivial correlated values of Vi , V2 for which the po- 
tential is non-dissipative (^(fc) =0 V fc) without being 
fully transparent. This set of values satisfies cosh(7r(7_|_) — 
cosh(7r(7_) 7^ 0, as can be seen from Eas. H26(l . Non-trivial 
solutions exist when g+,g- € C, yielding. 



1 = 71^1 - 4^1 - 4n2 nEl+. 



(31) 



Let us also notice from Eas. (|28|l that these solutions are 
also the borders where the physical scattering changes 
from the absorptive side to the emissive side or vice 
versa. We shall refer to these borders as inversion points 
(IP). Therefore, whenever we encounter an IP we can 
say ^(fc) = V fc without a fully transparent behaviour, 
and hence a non-dissipative scattering process for all en- 
ergies with a non-vanishing imaginary part of the po- 
tential. Let us note that from Ea. (|31|l the IP only ap- 
pears in the case of Scarf potential wells and only for 
1^2! < 1/4 — Vi. In Fig[7|the characteristic scattering 
probabilities are shown for a Scarf barrier and a Scarf 
well, and in Fig|Slthe maximum value of the physical ab- 
sorption is plotted versus IV2I for different values of Vi. 
When Vi is positive the absorption grows with the am- 
plitude of the imaginary part of the potential. On the 
other hand, for negative Vi a strikingly different pattern 
arises, with transparencies (T) and inversion points (IP) 
and the absorption does not increase monotonically with 
\V2\. The whole behaviour of the scattering can be clearly 
understood by building a two dimensional diagram IV2I 
vs. Vi (Figl^l, including physically permitted ranges, in- 
version lines and the points of transmission resonance. 
The complex Scarf potential shows two opposite faces to 
scattering, namely barrier and well, and a much richer 
structure in the latter case. 

After this detailed analysis of the peculiarities of the 
complex Scarf, that to our knowledge have not been re- 
ported before, let us continue with our work on connect- 
ing several potentials to model dissipative atomic chains. 
The procedure is the same as the one described for the 
Poschl- Teller potential hole in Section Only the por- 
tion of the potential included in the interval [—d, d] is 
taken (FiglS}. Making use of the asymptotic forms of 
the wave functions, and after some algebra, one finds the 
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FIG. 8; Maximum value of the absorption probability Amax 
vs. IV2I for different values of Vi. For each graph the up- 
per(lower) curve corresponds to a barrier(well). A dotted line 
is used in the forbidden values of | V2 1 , nevertheless the whole 
curve is shown for continuity. The allowed ranges of IV2I and 
the position of the IP can be compared with the values in 
Table El 



IV2I 




FIG. 9: Scattering diagram for the complex Scarf potential in 
terms of the potential amplitudes. The physically acceptable 
ranges for Vi, V2 correspond to the shaded zones. The curves 
are the inversion lines given by Eg. 1311 . The black points 
mark the correlated values of the amplitudes (Eqs.lPOJ) gen- 
erating a fully transparent behaviour. 



transmission matrix for the complex Scarf potential with 
a symmetrical cut-off, 



M(fc) 



■gi(v+2fed)^Y+li 



_^^~z{v+2kd)^l + SS 



(32) 

The matrix proves to be the same as for the asymptotic 
case (Ea. (|22() ). but with the extra phase 2fcc? in the di- 
agonal terms. It can be checked that the half-widths 
di^2 of both the real and imaginary parts coincide and 
that the decay of the imaginary part of the potential is 
always slower than that of the real part (Fig|Hl). This 
causes an increase in the minimum value of the cut-off 
distance do with regard to the Poschl- Teller case. For 
sensible values of the potential amplitudes we have found 
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k (x ^ units) k (x ' units) k (x ' units) 

FIG. 10: (Color online) Scattering probabilities for double Scarf potentials with parameters a = \[x~^ units), Vi = 2, V2 = 1 
for barriers and a = l{x~^ units), Vi = —4, V2 = 3.1 for wells. The solid lines were obtained from the analytical composition 
technique and the dashed lines correspond to the exact numerical integration of the Schrodinger equation. The insets show 
the potential profile (in arbitrary units) for each case: solid line for the real part and dashed for the imaginary part. The 
vertical lines limit the portion of the potential that the composition technique takes into account. For numerical integration the 
whole potential profile was considered. Notice that for low d the exact integration may lead to unphysical scattering because 
the conditions obtained for physical scattering are not valid for such d values. As d is increased an acceptable scattering 
is recovered and both methods start converging. Convergence is reached faster in the case of potential barriers where the 
analytical composition technique works impressively well even for very low values of d. 



that considering c?o — 3(ii/2 — 5.3/q! is enough in most 
cases. In fact, this minimum value can be relaxed in the 
case of potential barriers {Vi > 0), whereas for potential 
wells taking d below this value to apply the connection 
equations may sometimes distort the results. The cor- 
rect behaviour of the connection procedure for d > do 
can be observed in FigllUlwhere the scattering probabil- 
ities obtained upon integrating the Schrodinger equation 
numerically are compared with those given by the ana- 
lytical composition technique, for two potential barriers 
and two potential wells with different choices of the cut- 
off distance and with high values for IV2I. Notice how in 
the barrier example convergence between the methods is 
completely reached for d ~ 2(ii/2 whereas for potential 
wells a further step is needed because the convergence is 
much slower. 

Once the connection of potentials has been success- 
fully made, one should ask which are the ranges of the 
potential amplitudes that provide an acceptable physical 
scattering in this new framework. Analysis of this issue is 
very non-trivial and quite complex analytically, but also 
very important because it determines whether this model 



remains useful when considering atomic chains. First, in 
the case of two potentials we have observed that choosing 
each individual pair of amplitudes (Vi, V2) belonging to 
a physical range, and selecting the signs of the imaginary 
parts so that the physical faces of both potentials point 
in the same direction, then an acceptable scattering for 
the composite potential can always be recovered at least 
from one of the two possible orientations (both physical 
faces to the right or to the left). In other words, con- 
sidering that the incident particle always collides with 
the left side (it comes from —00) and therefore orientat- 
ing the individual physical faces to the left, then at least 
one of the sequences Vi{x) — Vii{x) or Vii{x) — Vi{x) 
gives an acceptable scattering for all energies. We have 
checked this assertion for a broad variety of Scarf couples. 
For a higher number of potentials the situation becomes 
more complex but a few pseudo-rules to obtain physi- 
cal scattering can be deduced. For an arbitrary chain 
we have found that in many cases the left scattering re- 
mains physical as long as: the left scattering of each in- 
dividual potential is physical and the left scattering of 
each couple of contiguous potentials is physical. This 
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recipe seems completely true when composing potential 
barriers only, whereas when wells are included it fails in 
some situations, especially when several contiguous wells 
are surrounded by barriers. Although at the beginning 
it may appear almost random to recover a physical scat- 
tering from a large composition of Scarfs, following the 
given advices it turns out to be more systematic. 

Let us remember that the composition procedure, 



apart from being a powerful tool for numerical calcu- 
lations also provides analytical expressions for the scat- 
tering probabilities, which of course adopt cumbersome 
forms for a large number of potentials but are useful for 
obtaining simple expansions for certain energetic regimes. 
Just as an example, the transmitivity and left reflectivity 
for the double Scarf read. 



2Scarf{ - ^2^2 + (1 + SiSi)(l + S2S2) - 25152^1 + ^iSl ^1 + S2S2 COs((pi + (^2 + 2A:(di + ^2)) 

/^x _ + ^2^2) +^2(1 + ■Sl^l) ~ 2S1S2VI + Sl5l\/1 + 52^2 COs(y?i + (^2 + 2fc(rfi + ^2)) ,„^x 

2Soar/l " ^2^2 ^ _^ SiSi)(l + S2S2) - 2siS2Vl + SlSi ^1 + S2S2 COs(^i + (^2 + 2/s(rfi + ^2)) 



making use of the previously defined terms s,s and ip. 
One important feature of the formulae for the composite 
scattering probabilities is the fact that they analytically 
account for the fully transparent behaviour of the whole 
structure as long as there is resonant forward scatter- 
ing of the individual potential units, as can be seen from 
the latter equations and also from Ea. H18|l for the dou- 
ble Poschl- Teller. Another curious feature arises when 
composing different potentials whose amplitudes describe 
an IP. In this case the whole structure remains non- 
dissipative fFig lll|) . Moreover, the complex Scarfs at 
these points behave completely as real potentials, pro- 
viding an acceptable scattering for the composition that 
is independent of the incidence direction for any sequence 
of the individual Scarfs. 

Considering larger Scarf chains with small imaginary 
parts of the potentials, in the case of a periodic array 
we observe that the absorption term remains flat over a 
wide range of forbidden bands and oscillates inside the 
permitted ones. The variations in the absorption are 
entirely balanced by the reflectivity while the transmi- 
tivity is surprisingly not affected by the presence of a 
small complex potential ('Fig ll2|l . This behaviour con- 
trasts strongly with the complex delta potentials peri- 
odic chain where the absorption was completely different 
and it was the reflectivity that was little affected by the 
dissipation. For an aperiodic sequence the situation is 
quite different, as expected. In Fig^Ja type of molecu- 
lar aggregate is modelled with complex Scarfs. It exhibits 
a peaky absorption spectrum and a strongly oscillating 
transmitivity with sharp resonances. In Fig ll4l a sym- 
metric atomic cluster has been considered in which the 
dissipation only occurs at both ends. Different transmis- 
sion and absorption configurations can be obtained by 
building different structures. This shows the usefulness 
and versatility of this model for being able to account for 
a variety of possible experimental observations. 

As a final exercise, two examples are included in Fig llSI 
showing the evolution of the transmitivity of two different 




k (x ' units) 



FIG. 11: (Color online) Composition of Scarfs with differ- 
ent pair amplitudes describing IPs. Sequence with a — 
l{x-^ units) and {Vi,V2): (-4,%/T3), (-4.5,^15), (-5,717), 
{—6,y/21), (— 7, 2vT^). AH the cut-off distances have been 
taken equal to d = 6(2: units). The inset shows the potential 
profile with solid(dashed) line for the real(imaginary) part in 
arbitrary units. Notice that the composition remains non- 
dissipative. 



Scarf compositions as a function of the imaginary part of 
the potential. The transmission patterns are plotted for 
different values of the parameter e = |V2|/|Vi|, which 
measures the strength of the imaginary part. The trans- 
mission efhciency is evaluated using an averaged trans- 
mitivity Tint = J^^ T{k) dk/{k2 — ki) corresponding to 
the area enclosed by T(k) per energy unit in a character- 
istic energy range, namely the zone where T evolves un- 
til it becomes saturated. The imaginary potential tends 
to smooth the transmission pattern in the first example 
(corresponding to a sequence of barriers) , causing a slight 
decrease in the averaged transmitivity for low e, although 
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k (x ' units) 



FIG. 12; (Color online) Periodic chain of 20 complex Scarfs 
with amplitudes a = l{x~^ units), Vi = 2, V2 ~ 0.02 (solid 
lines) and V2 = 0.05 (dashed lines) and cut-off distance d = 
6 (a; units). The transmitivity is not changed much by the 
small imaginary part of the potentials. 




k (x ' units) 



FIG. 13: (Color online) Scattering probabilities for 
a 10-Scarf structure with parameters a — l{x~^ 
units) and {Vi,V2): (1,0), (-0.5,-0.005), (-1.3,0.013), 
(1.8,0.01), (-1.3,0.013), (4,0.04), (-2.4,-0.024), (3.5,0.04), 
(-3.1,0.031), (-1.5,0.04). The cut-off distance is d = 6{x 
units) equal for all of them. The inset shows the real part of 
the potential profile (arbitrary units). 



it is finally improved. For a double well a different effect 
occurs. Tint is always enhanced with increasing e un- 
til it reaches a maximum, after which the transmitivity 
falls with the imaginary potential. With these results in 
mind one could speculate that these models might also 
be useful to describe the phase-breaking inelastic scat- 
tering processes in atomic chains. Especially when one 
treats disordered arrays where the break of the localiza- 
tion regime could arise as a result of the loss of coherence 
due to inelastic collisions. 




k (x ' units) 



FIG. 14: (Color online) Scattering probabilities for a 10- 
Scarf structure with parameters a = l{x~^ units) and 
iVi,V2): (1,0.04), (1,0.05), (-2,0), (-1.5,0), (-1,0), 
(-1,0), (-1.5,0), (-2,0), (1,0.1), (1,0.1). The cut-off dis- 
tance is d = 6(x units) equal for all of them. The inset shows 
the real part of the potential profile (arbitrary units). 

VII. CONCLUDING REMARKS 

In this work we have used the transfer matrix for- 
malism in the form of the multiple collision technique 
to model dissipative scattering processes by using com- 
plex potentials in various models of atomic chains. The 
absorption probability has been rigorously included to 
recover unitarity for the non-hermitian hamiltonians. 

New exact analytical expressions are given for the scat- 
tering amplitudes of an arbitrary chain of delta poten- 
tials. The absorption effects arising by promoting the 
delta couplings to the complex domain have been shown, 
revealing the flexibility of this simple model to account 
for very different physical schemes. 

Handy expressions for the transmission matrices of the 
Poschl- Teller and the complex Scarf potential have been 
constructed in their asymptotic as well as their cut-off 
versions for the first time. These latter matrices have 
made feasible, via the composition technique, the assem- 
bly of an arbitrary number of potentials to build atomic 
one-dimensional wires that can incorporate absorptive 
processes. Different absorption configurations are pre- 
sented in several examples that show the versatility of 
the model to account for a variety of possible experimen- 
tal observations. The procedure developed is not only 
useful for numerical calculations but also provides ana- 
lytical formulae for the composite scattering probabili- 
ties, which have not been obtained by other methods, 
and whose significance has been checked by numerical 
integrations of the Schrodinger equation. 

From a more mathematical view, a complete and rig- 
orous analysis of the scattering properties of the com- 
plex Scarf potential has been carried out. The ranges of 
physical transmission have been obtained and a group of 
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FIG. 15: (Color online) Transmitivity for different strengths of the imaginary part of the potentials in (a) a chain of 6-Scarf 
barriers Vi — 4,2,2.5,3,3.5,4 and a = l{x~^ units), (b) a double Scarf well with Vi = —1,-1 and a — 2{x~^ units). The 
cut-off distances in all cases are equal to d = 6(2; units). Each imaginary amplitude reads V2 = except the first barrier 

in example (a) which is maintained real. The insets show the evolution of Tint vs the strength of the imaginary amplitudes. 



novel features have arisen such as the presence of perfect 
transparencies and inversion points. 



APPENDIX A: SOLUTIONS OF THE 
POSCHL-TELLER POTENTIAL 



Apart from being able to include dissipation in the sys- 
tems in a tractable way, the tools and methods provided 
may have direct applicability for considering molecular 
aggregates |2^ and other structures with explicit poten- 
tial profiles and also to build a new kind of continuous 
disordered models. 



Future work seems promising since we are already in 
progress on the assembly of the one-dimensional struc- 
tures to quantum dots in order to analyze the effect of 
dissipation on the conductance. Our aim is also to treat 
long disordered wires within this framework, and in par- 
ticular to study the applicability of these potentials to 
account for the loss of coherence due to inelastic scatter- 
ing process in the electronic localization regime. 



The elementary positive energy solutions of the 
Schrodinger equation with potential H14|l read, 



e{x) — cosh'^(Qa;)F ^a, 2 ' ~ sinh^(aa;) 

o{x) =sinh(ax) cosh^(Q;x)x 
/I 13 



(Ala) 



(Alb) 



where a — ^ + i^, b — ^ — and F(a, b, c: z) is the 

2 2a ^ 2 2a \ 7 : : / 

Hypergeometric function. And their asymptotic forms 
can be written as, 



(x) sign(a;) (^ne 



— iO —ik\x\ I _ iO ik\x\ 



where 



V^r(z|)2- 



ne 



V^r(z^)2-# 



(A2a) 
(A2b) 



(A3a) 
(A3b) 



2a J 
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1. Scattering states 

The elementary positive energy solutions of the 
Schrodinger equation with the complex Scarf potential 
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(O are, 



2. Physical transmission 



/ k k I i 

F \ c + i — , c ~ i — , 1 — b + c] — + — sinn(aa;) 
\ a a 2 2 



(Bla) 



- + ^sinh(aa;)j x 

/ k k 1 j-wN 

F \ b + i—,b — + b — c; — \ — sinn(ax 

\ a a 2 2 

(Bib) 

F{a, b, c; z) being the hypergeometric function and using 
the 6, c definitions in H25|l . The asymptotic hmit x — > ±00 
yields, 

mix) ^ 2^= [sign(x) • (u;ie-''=l^l + ^26**^1^1) 

(B2a) 

U2{x) -^-2" [-sign(a;) -i] 5 (gie^'^^l^l + gse'*^'''!) 

(B2b) 

where 

^"1 = — ..L.. ?^ (B3a) 



W2 = 



T{c-^>t)T{l-b-^^t) 
2-2»4r (1 - 6 + c) r (2i|) e"'s"(^)f ■ 



T{c + ^^t)T{l-b + ^^t) 
qi ^wi{b^ c) 
q2 ^W2{b^c). 



(B3b) 

(B3c) 
(B3d) 



The condition for physical transmission is 
cosh(7r(7+) cosh(7r(7_) > 0. V2 can be considered 
positive with no loss of generality since its change in sign 
(which is equivalent to changing the side of incidence) 
does not affect the transmission. With the definitions 
X = I Vi I — V2 and 1" = I Vi I + V2 , the study can be easily 
carried out. Considering Vi > the inequality translates 
into the permitted regions 



V Y 



X>0 

-2n{2n + I) < X < -2n(2n - 1) n G Z+ 



(B4) 

which is clearly a sequence of allowed vertical fringes in 
the negative X quadrant and the whole positive X quad- 
rant. This pattern will be same but rotated 7r/4 clockwise 
when the change of variables is undone. More specifically, 
in terms of the potential amplitudes the allowed intervals 
can be written as 



\V2\ e[o,i/i]u 

[2n(2n- 1) + Fi,2n(2n+ 1) + Fi] n G Z+. 

(B5) 



In the case of Vi < a careful analysis leads to the 
following cumbersome allowed sets 



X 2n{2n - 1) < F < 2n(2n + 1) n = 1, 2, 3, . . . 

f{2n(2n- I) < X < 2n(2n + 1)} f] {2m{2m -1) <Y < 2m(2m+ 1)} n,m = 1,2, . . . 
j{2j(2j + 1)<X < 2j(2j + 3) + 2}n {2fc(2fc+ 1) <Y < 2fc(2fc + 3) + 2} k,j = 0,1,2,. 



X >0 



(B6) 
(B7) 



which are a set of allowed horizontal fringes in the nega- 
tive X quadrant and a chessboard like structure for posi- 
tive X. Undoing the change of variables will mean a 7r/4 
clockwise rotation followed by a reflection around the ver- 
tical axis of this pattern to recover the negative axis of 



Vi. Solving these inequalities in terms of the potential 
amplitudes gives rise to the following set of inequalities, 
each one assigning a certain allowed interval for | V2 1 when 
fulfilled. 



for n e Z+, 



\Vi\< n(2n- 1) 
n(2n- 1) <\Vi\< n(2n+ 1) 



[2n(2n- 1) - \Vi\ , 2n{2n + 1) - \Vi\] 
[\Vil2n{2n+l)~\Vi\] 



(B8) 
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for m,n £ Z"*" m > n, 

m{2m- 1) + n(2n- 1) <\Vi\< m(2m - 1) + n(2n+ 1) ^ [2m(2m - 1) - \Vi\,-2n{2n- 1) + |yi|] 

m{2m - 1) + n(2n+ 1) < \Vi\ < m{2m + 1) + n{2n - 1) ^ [- 2n(2n+ 1) + |yi|,-2n(2n- 1) + l^il] (B9) 

m{2m + 1) + n{2n — \) <\Vi\ < m(2m + 1) + n{2n + 1) [- 2n(2n+ 1) + \Vil2m{2m + I) - \Vi\] 

and for j, /c = 0, 1, 2, . . . fc > j, 

fc(2fc + l)+j(2j + l) < l^il <fc(2fc + l)+j(2j + 3) + l^ [2k{2k+\)-\Vil-2j{2j + \) + \Vi\] 
k{2k + 1) + j(2j + 3) + 1 < l^il < k(2k + 3) + 1 + j(2j + 1) ^ [ - 2 - 2j(2j + 3) + |Vi|, -2j(2j + 1) + l^il] 
k{2k + 3) + 1 + ]{2j + 1) < l^il < k{2k + 3) + ]{2j + 3) + 2^[-2 - 2]{2j + i) + \Vi\,2 + 2k{2k + 3) - |Vi|] . 

(BIO) 

I 

In the particular cases m — n for Eas. (|B9|) and k — j for must be considered. The total physical range for IV2I 
Eas. ljB10|) only the positive part of the allowed intervals comes from the union of the different permitted intervals. 
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